Magnetothermal transport in the spin- 1/2 chains of copper pyrazine dinitrate 
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We present experiments on the thermal transport in the spin- 1/2 chain compound copper pyrazine 
dinitrate Cu(C4H4N2)(NOa)2. The heat conductivity shows a surprisingly strong dependence on 
the applied magnetic field B, characterized at low temperatures by two main features. The first one 
appearing at low B is a characteristic dip located at [ibB ~ fcsT, that may arise from Umklapp scat- 
tering. The second one is a plateau-like feature in the quantum critical regime, (j,b\B — B c \ < UbT, 
where B c is the saturation field at T — 0. The latter feature clearly points towards a momentum 
and field independent mean free path of the spin excitations, contrary to theoretical expectations. 

PACS numbers: 75.40.Gb 66.70.+f, 75.47.-m 



In the last decade, considerable progress has been 
achieved in theoretical studies of thermal transport in 
one-dimensional (ID) quantum spin systems (for a recent 
review see Ref. 1). One of the most important model sys- 
tems is the Heisenberg spin 5=1/2 chain with isotropic 
antiferromagnetic interactions, described by 

N N 

H = jJ2S 1 S 1+1 -9VbBY / SI (1) 
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where J is the intrachain nearest-neighbor exchange and 
g/j,B the magnetic moment. Recently, a number of non- 
trivial effects were predicted for the spin thermal con- 
ductivity (k s ) of this system in external magnetic fields 
£ 2,3,4,5,6,7 of particular interest is the behavior close to 
the saturation field B c = 2J/gfiB, which defines a quan- 
tum critical point. 8 Experimental information on k s (B) 
is, however, missing because most of the known realiza- 
tions of the Hamiltonian (1) have large values of the in- 
trachain exchange constant J/ks ~ 100 — 1000 K, that 
severely limits the region of the phase diagram accessi- 
ble using standard laboratory equipment. The existing 
results are limited to studies of the phonon thermal con- 
ductivity n ph {B, T) in CuGe0 3 and Yb 4 As 3 , 9 ' 10 ' 1M2 and 
do not address the behavior of k s (B, T). 

Copper pyrazine dinitrate Cu(C4H 4 N2)(N03)2 
(CuPzN) appears to be an ideal compound for the 
thermal conductivity experiments in magnetic field, as 
J/k B = 10.3 K and therefore B c = 15.0 T, 13,14 which 
is easily accessible experimentally. CuPzN has an or- 
thorhombic structure with lattice constants a = 6.712 A, 
b = 5.142 A, and c = 11.73 A at room temperature. 
The chains of Cu 2+ spins 5=1/2 run along the a axis. 
Inelastic neutron scattering, magnetization, and specific 
heat measurements have confirmed that CuPzN is very 
well described by the model of Eq. (l).i3,U,i6 The 1D 
nature of the spin interaction is reflected by a very low 
ordering temperature, Tjv = 0.107 K, and therefore the 
ratio of interchain (J') to intrachain (J) couplings is 
estimated to be tiny, \J'/J\ « 4.4 x 10~ 3 . 17 



In this paper, we report measurements of the thermal 
conductivity n(B, T) of CuPzN in the temperature re- 
gion between 0.37 and 10 K and in magnetic fields up 
to 17 T. The crystals of CuPzN were grown from water 
solution of pyrazine and Cu nitrate via slow evaporation. 
The crystals have right-prism shapes with the a axis of 
length 10 mm directed along the height of the prism. 
The dimensions perpendicular to the a axis are typi- 
cally of the order of 0.4x0.7 mm 2 . Thermal conductivity 
was measured by a standard steady-state heat-flow tech- 
nique, where the temperature difference was produced 
by a heater attached to one end of the sample and mon- 
itored by a matched pair of RuC>2 thermometers. The 
temperature difference between the thermometers was of 
the order of 1% of the mean temperature T. In one sam- 
ple the heat flux was oriented along the a axis to measure 
the thermal conductivity parallel to the chains (/«"). On 
another sample, the thermal transport perpendicular to 
the chains (k ± ) was measured. Magnetic fields were ori- 
ented parallel to the a axis. 

In Fig. 1, we show the temperature dependence of k" 
and k 1 - in constant magnetic fields. The important ob- 
servation is that B strongly influences k", but no signif- 
icant changes of k with field are observed. In CuPzN, 
the magnetic interaction between the spin chains is ex- 
tremely weak. Therefore one expects that the heat cur- 
rent perpendicular to the chains is of purely phononic 
origin, k 1 - ea n^ h , consistent with the observed absence 
of a significant field dependence. Along the chain di- 
rection, k" = Kp[j + K m , one can separate the field- 
independent phononic contribution, Kp\(T), unrelated to 
the presence of the spin chains, from a B dependent mag- 
netic contribution, K m {B,T). Three terms contribute to 
K rn = Aftp h + k s + K sp . First, the scattering by spin 
excitations affects the lifetime of phonons and gives rise 
to a decrease AKp h of the phononic thermal conductiv- 
ity. Second, there is a spin contribution k s to the ther- 
mal transport. Finally, a spin-phonon cross term K sp de- 
scribes the 'drag' of spin heat currents by phononic heat 
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FIG. 1: The thermal conductivity of CuPzN parallel and 
perpendicular to the spin chains as a function of temperature 
in several constant magnetic fields. The solid circles corre- 
spond to the calculated zero-field spin thermal conductivity 
along the chains (see text). 



currents (and vice versa). For constant T, the difference 
between the total measured ac" (B) and its zero-field value 
At" (B = 0) allows one to extract the field dependence of 
the magnetic heat conductivity n m {B, T) — ac to (0, T). 

In Fig. 2, we plot [k)\(B) - k)\(Q)]/T w [n m {B) - 
K m (0)]/T at several constant temperatures. Depending 
on T, one can distinguish two types of behavior. The 
crossover between these two types occurs at about 2.5 - 
3 K. At low temperatures, T <C J/ks, K m (B) has two 
main features: a decrease with a minimum at low fields 
and another decrease at higher fields. On top of the high- 
field decrease, there is a plateau-like feature in the vicin- 
ity of the critical field B c . At high temperatures, K m (B) 
has a simpler shape with a single minimum close to B c /2. 
The circles in Fig. 3 show the positions of the K m (B) min- 
ima. In the same figure the triangles, determined as the 
inflection points of n m {B) curves, characterize the size of 
the plateau- like regions in the vicinity of B c . 

For the interpretation of these data, we first note 
that a spin gap = g^isiB — B c ) opens when B ex- 
ceeds B c . 18 Therefore, for gfi B {B - B c ) > fc B T, both 
heat transport by spin excitations and phonon scatter- 
ing by spin excitations should vanish. In our exper- 
iment, K m (B) decreases with increasing B > B c for 
low T < 2.5 K. Thus, the observed field dependence of 
the thermal conductivity arises dominantly from ther- 
mal transport in the spin system and not from AAtJ| h . 
The overall decrease of At m towards lower T, see Fig. 1, 
suggests that impurity scattering dominates spin-phonon 
scattering. Therefore, spin-phonon drag terms should 
not be important. This allows us to associate n m with 
the spin thermal conductivity k s . We model k s using 
a combination of a relaxation time approximation and 
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FIG. 3: A (B,T) diagram of the characteristic features of the 
thermal conductivity of CuPzN. The circles correspond to the 
minima of K m {B). The triangles correspond to the inflection 
points of n m (B) in the vicinity of B c . The solid lines are 
the positions of the K m {B) inflection points calculated for a 
constant mean free path (see text). 



a mean-field theory (MFT) for the spin chain relying 
on mapping Eq. (I) onto a system of interacting spin- 
less fermions via the Jordan- Wigner transformation. 19 
The total weight of the frequency-dependent spin ther- 
mal conductivity, / k(uj)cLj, in Heisenberg chains calcu- 
lated within the MFT approach has been shown to be 
in reasonable agreement with exact results of the Bcthe 
ansatz calculations. 6 ' 20 The fermions occupy a cosine 
band which, at B = 0, is half-filled with the chemical po- 
tential at kp = ±7r/2a. An external magnetic field shifts 
the chemical potential and changes the band width. The 
dispersion of the fermions is 

e k = - J(l + 20) cos(ka) - g/j, B B + 2Jm, (2) 

where k is the wave vector and a is the distance between 
neighboring spins. The parameter and the average 
local magnetization m are determined self-consistently 
from O = f H I<X cos(ka)f k dk and m = - | + f J„ /a fkdk 
where f k = (cxp(£fe/fcsr) + I) _1 is the Fermi distribution 
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function. The spin thermal conductivity is given by 
Na T 

K s = L2, 



T f ^ dfk n -l 2 i) 

Ln = J o ~^e k V k T k dk, 



(3) 
(4) 



TV is the number of spins per unit volume, v k = dek/dk 
is the velocity and r k is the relaxation time related to the 
mean free path l k as r k — l k /v k - 

As discussed above, scattering by defects is expected 
to be the most important source of extrinsic scatter- 
ing at low T. Assuming for the moment a mean free 
path, lk(B, T) = l(T), which is independent of both mag- 
netic field and momentum, we can calculate k s ^ _1 from 
Eq. (3) as a function of B and T without free param- 
eters. For our calculations, we used the experimental 
values of the g-factor parallel to the a-axis g = 2.05 from 
Rcf. 21 and J/k B = 10.3 K from Rcf. 14. The result- 
ing [k s (B) — k s (Q)]/IT is shown in Fig. 4(a) for several 
T. The calculations obviously give correct qualitative 
account for the experimentally observed behavior at in- 
termediate and high fields for all T. It reproduces the 
crossover from the low-T to the high-T behavior around 
3K and predicts the low-T plateau-like features at B c 
with the correct width. In Fig. 3, the calculated inflec- 
tion points of n s {B) (solid lines) are in good agreement 
with the experiment. Close to the saturation field, i.e. at 
the quantum critical point, one finds the following scaling 
relation 
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where f(x) is a dimensionless scaling function indepen- 
dent of the precise dispersion of the spin excitations. 

The comparison of the theoretical curves for various as- 
sumptions on the momentum dependence of the scatter- 
ing rate, Fig. 4(b), with our experimental results, Fig. 2, 
clearly supports a mean free path which depends on T 
but not on momentum and magnetic field at least close 
to the quantum critical point. Note that close to B c , the 
mapping of the spin chain to free fermions becomes ex- 
act and one may therefore hope that a simple Boltzmann 
description becomes asymptotically exact (if strong local- 
ization by disorder can be neglected). Here one has to 
emphasize, that the observation of a constant mean free 
path is highly surprising in a one-dimensional system. In 
the limit of weak disorder (small potential strength com- 
pared to T or \B — B c \) the golden-rule scattering rate 
1/Tfc is proportional to the density of states, 1/rfe oc l/v k , 
where v k is the velocity of the spin excitations, and there- 
fore one gets l k cx v 2 . Luttinger liquid corrections to 



this formula are expected to vanish close to the quan- 
tum phase transition when the density of spin excitations 
is small. In the opposite limit of strong impurities, the 
spin excitations have to tunnel through the potential and 
the conductivity is proportional to the transmission rate, 



ffc oc v. 



implying that I oc T k oc v 2 ., as in the limit of 



weak disorder. However, the heat conductivity, calcu- 
lated assuming l k oc v 2 , is featureless at B « B c , in clear 
disagreement with the experiment, see Fig. 4(b). 

What are possible mechanisms which can lead to a con- 
stant mean free path? In the absence of inelastic scat- 
tering, interference effects lead to localization of the spin 
excitations but neither the observed T dependence of k, 
see Fig. 1, nor the constant mean free path point towards 
the importance of localization effects. Nevertheless, rare 
defects may cut the spin chains into separate pieces of 
finite length. In such a situation, the effective mean free 
path is determined either by scattering from one spin 
chain to the next or - more likely - by coupling heat cur- 
rents in and out of these pieces of spin chains by lattice 
vibrations. The first scenario naturally leads to a con- 
stant mean free path, but a theoretical prediction for the 
more likely second scenario is presently lacking. 

Recently, several theoretical papers considered magne- 
tothcrmal corrections (MTC) to the spin thermal con- 
ductivity, which should appear in a magnetic field due to 
a coupling of the heat and magnetization currents. 2,4,5,6 
In this case, n s is given, instead of Eq. (3), by 6,20 



Na 



(6) 



The second term in parentheses represents the MTC. 
However, in Ref. 3 it was argued that MTC are absent in 
macroscopic samples of real materials where the conser- 
vation of the total magnetization parallel to B is broken 
by spin-orbit coupling, prohibiting a piling up of magne- 
tization. In Fig. 4(b), we show k s (B) calculated within 
the relaxation time approximation for one temperature 
without (solid line) and with (dashed line) MTC. Ac- 
cording to these calculations, MTC should completely 
destroy the plateau-like feature at B Cl obviously at vari- 
ance with the experiment. Thus, our experiment provides 
strong evidence against the existence of magnetothermal 
corrections to the spin thermal conductivity in CuPzN. 

The assumption of a constant mean free path agrees 
well with the experiment at intermediate and high fields 
but predicts a field-independent k s at low fields for 
T <C J/ks- The experiment, however, shows a low-field 
minimum of n s (B) located approximately at [i B B ~ k B T 
for low T, see Fig. 3. This suggests strongly that the ori- 
gin of this anomalous behavior is related to the physics 
of the half-filled band of Jordan- Wigncr fermions. Only 
excitations in a window of width k B T around the B = 
Fermi surface are able to relax their momentum to the 
lattice by Umklapp scattering. An explanation of the 
minimum might be possible along the following lines. 
The spin excitations which dominate the thermal trans- 
port have a typical energy of order k B T, those located 
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FIG. 4: (a) Spin thermal conductivity as a function of B 
for several T, calculated for CuPzN within a relaxation time 
approximation assuming a B- and momentum-independent 
mean free path I. (b) Spin thermal conductivity at T = 
0.66 K, calculated under various assumptions for the mean 
free path (see legend), and taking into account magnetother- 
mal corrections (dashed line) . Only the model with constant 
I is consistent with the data. 



directly at the Fermi surface with energy do not con- 
tribute. Therefore, Umklapp scattering is most effective, 
if the excitations with energy fc^T (rather than 0) have 
a momentum close to ±7r/(2a), where Umklapp scatter- 
ing is strongest, i.e. for [IbB ~ ksT. Indeed, if one 
assumes a scattering rate with a peak at k = ±7r/(2a), 
the relaxation time approximation yields a minimum in 
k(B) at the right position (not shown). An exponentially 
strong maximum in k(B) has been theoretically predicted 
for clean spin chains coupled to phonons in Ref. 3 as a 
consequence of the existence of certain approximate con- 
servation laws. As our system is disorder dominated, we 
do not believe that this mechanism is directly applicable 
in the present situation. A minimum in k s (B) was also 



predicted for a classical ID Heisenberg model in Ref. 7 
but with a different T dependence. A full theory of the 
low-field minimum should include the interplay of Umk- 
lapp, impurity and phonon scattering which is currently 
under investigation. 

According to Eqs. (3, 5), at low T, k s (B c ,T) ~ 
k s (0,T)/2, see Fig. 4(b). Using this, we estimate the 
absolute values of the zero-field spin contribution as 
k s (0,T) w 2[k(0,T) - k(B c , T)], which is shown in 
Fig. 1. It is notable that n s oc T 2 implying a linear 
increase of the mean free path with T as I = AT, with 
A = 1.0 x KT 6 m/K for T < 1.5 K. This observation is in 
agreement with theoretical predictions for weakly disor- 
dered spin chains, where the linear T dependence arises 
from the renormalization of the impurity potentials by 
Luttinger liquid corrections, see e.g. Refs. 22,23. 

In summary, we have experimentally established the 
magnetic contribution n m to the thermal conductivity 
of the 5 = 1/2 chain compound copper pyrazine dini- 
trate. At low temperatures, the field dependence of n m 
is characterized by two features, one at high fields in the 
vicinity of B c and the other at low fields. The low- field 
feature remains to be explained. The high-field feature is 
associated to the spin thermal transport k s with a mean 
free path surprisingly weakly dependent on both field and 
momentum. No magnetothermal corrections to k s have 
been identified in our experiment. 
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